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Abstract We report on the realization of a double- 
well potential for Rubidium-87 Bose-Einstein conden- 
sates. The experimental setup allows the investigation of 
two different dynamical phenomena known for this sys- 
tem - Josephson oscillations and self-trapping. We give 
a detailed discussion of the experimental setup and the 
methods used for calibrating the relevant parameters. 
We compare our experimental findings with the predic- 
tions of an extended two-mode model and find quanti- 
tative agreement. 

PACS: 74.50+r, 03.75.Lm, 05.45.-a 



1 Introduction 

One of the most prominent features of quantum me- 
chanics is the tunnelling of massive particles through 
classically forbidden regions. Although tunnelling is a 
purely quantum mechanical effect, it can be observed 
on a macroscopic scale if the system can be described 
by two weakly linked (i.e. having a small spatial over- 
lap) macroscopic wave functions with global phase co- 
herence. A fundamental physical phenomenon based on 
macroscopic tunnelling is the Josephson effect predicted 
by Brian D. Josephson in 1962 1 . The first observation 
of this effect has been realized with two superconductors 
separated by a thin insulating barrier and was reported 
by Anderson et al. 2 just one year after its theoretical 
prediction. Since then many experiments with the elec- 
tronic Josephson junction system have been performed. 
The Josephson effect has also found its way to techno- 
logical applications such as superconducting quantum 
interference device (SQUIDS) which allow to measure 
weak magnetic fields with very high precision. 

The Josephson effect for neutral superfluids has al- 
ready been observed with liquid ^He 3 and ^He 4 ex- 
hibiting the typical current-phase relation. It is charac- 
terized by an alternating current flowing through the 



central tunnelling barrier if a constant energy differ- 
ence between both sides is applied. In the context of 
dilute quantum gases such as Bose-Einstein condensates 
(BECs) Josephson junction arrays have been demon- 
strated ini and only recently a single weak link was re- 
alized and the predicted tunnelling dynamics has been 
observed 0. The main subject of this paper is the discus- 
sion of this experimental setup. We will give a detailed 
description of the experimental implementation and the 
necessary thorough calibration of system's parameter. 
We are also going to discuss the comparison of the ob- 
tained data with an extended two-mode model recently 
developed by D. Ananikian and T. Bergeman [Zj. 



2 Basic setup 

There are many different methods to produce a double- 
well potential for Bose-Einstein condensates. The first 
realized double-well potential was obtained already in 
the early days of BEG using a harmonic confinement 
created by a magnetic trap and a focussed blue detuned 
laser beam generating the tunnelling barrier. This setup 
with a well separation of typically 50/im was used for 
the first interference measurements with BECs 8 . Many 
other attempts have been made to realize double-well 
potentials |9lll()mi but so far only the splitting of con- 
densates into two independent parts was possible. Re- 
cently we succeeded to observe Josephson dynamics in a 
double- well potential realized with optical dipole poten- 
tials 6 . The basic idea is the combination of a three di- 
mensional harmonic confinement and a one dimensional 
periodic potential with a large lattice spacing of 5 /im 
(see fig.nj. 

The ^^Rb BEC in our experiment is prepared in a 
crossed beam dipole trap consisting of two orthogonal 
far red detuned Nd-Yag laser beams (see figQ). A fre- 
quency difference of lOMHz between the two beams is 
introduced in order to avoid uncontrolled interferences 
between both beams. The dipole trap beam is radially 
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Fig. 1 Realization of a single Josephson junction using opti- 
cal dipole potentials, a) The basic setup is given by a crossed 
beam dipole trap combined with a periodic light shift po- 
tential implemented with two beams intersecting under an 
angle of 9°. b) On the scale of the three dimensional har- 
monic confinement the periodic potential is hardly visible, c) 
In the center of the combined potential a clean double-well 
potential is realized. 



symmetric and has a waist of 60(5)/im. It provides con- 
finement in the direction of gravity (in the following 
discussion z-direction) and in y-direction. The crossed 
dipole beam is radially asymmetric and has a radius at 
the intersection of 140(5)/im in z-direction and 70(5)/im 
in x-direction. This asymmetry allows the adjustment 
of the harmonic confinement along the x-direction, i.e. 
the direction of the double- well potential, without signif- 
icantly changing the trap frequencies in the other direc- 
tions. With a maximum power of approximately 500mW 
in the dipole trap beam and 800mW in the crossed dipole 
beam, we can achieve a maximum three-dimensional con- 
finement of bJx.max ^ 27r X 120(6)Hz, UJy^rnax ~ 27r X 

170(8)Hz and ijJz.max ^ 27r x 180(8)Hz and a maximum 
trap depth of approximately 5/iK. The typical sponta- 
neous emission rate is estimated to be O.OlHz and thus 
is negligible for all our experiments. 

The periodic potential is realized by a pair of laser 
beams with parallel linear polarization crossing at a rela- 
tive angle of Qf=9(l)° as depicted in fig.^a). The beams 
are aligned symmetrically with respect to the crossed 
dipole beam, such that the resulting potential is modu- 
lated only in x-direction. The laser wavelength is chosen 
to be 811nm. The resulting potential is shown in fig. 
^b,c) and is given by 

Vdw = ^rnwlix - Axf + Vb cos^ ^^ j (1) 

where Ax is the relative position of the two potentials. 
For Ax = a symmetric double-well potential is ob- 
tained. 



The intensities of all laser beams and the relative 
distance between the center of the crossed dipole beam 
and a maximum of the periodic potential are very criti- 
cal parameters thus an active stabilization is inevitable. 
The light intensities are stabilized to better than 10~^ 
which allows the generation of small condensates with 
small shot-to-shot fluctuations of the final atom num- 
bers, i.e. 1100(300) atoms. The control of the position 
of the crossed dipole beam with respect to the standing 
light wave is crucial for the investigation of the dynam- 
ics of bosonic Josephson junctions. The relative shift de- 
termines the resulting shape of the effective double- well 
and its symmetry. In order to stabilize the position of 
the periodic potential, the relative phase of the two lat- 
tice laser beams is controlled with two acousto-optical 
modulators. The absolute position of the periodic po- 
tential maxima is deduced by coupling out 2% of the 
beam intensities and bringing them to interference un- 
der an angle of 6(1)°. This leads to an interference pat- 
tern with a periodicity of about 8/im. A precision air 
slit with a width of 5(1) fivn and a height of 3mm is 
placed in the overlap region and is oriented parallel to 
the interference fringes. The light intensity transmitted 
through the slit, which directly depends on the relative 
phase between the two beams, is monitored with a pho- 
todiode. The signal of the photodiode is then locked to a 
given value using a standard proportional-integral loop 
amplifier. The achieved position stability is deduced by 
analyzing the position of the split-BEC images and gives 
an upper bound for the stability. We obtain a standard 
deviation of approximately lOOnm (limited by the pixel 
size of the camera system), which corresponds to a phase 
stability of less than 7r/50. 

In order to initially prepare a given population imbal- 
ance, we use a controlled shift Ax of the crossed dipole 
trap beam which results in an asymmetric double-well 
potential, where more atoms are accumulated in the 
lower well. This shift is experimentally implemented with 
a piezo actuated mirror mount. The achievable popula- 
tion imbalances and the associated uncertainties are de- 
picted in fig.EJ For Z\x=0 the ramping up of the periodic 
potential into the condensate prepared in the harmonic 
potential leads to a symmetric population of both wells 
with vanishing population difference z = {Ni—Nr)/{Ni-\- 
Nr) = {Ni^r number of atoms in the left and right well, 
respectively). A shift of only Ax = 500nm already intro- 
duces a population imbalance of Zc = 0.39. As we will 
show in the last section this corresponds to the critical 
imbalance for our experimental parameters which sepa- 
rates the regime of Josephson oscillation dynamics, i.e. 
z < Zc (indicated with gray shading in fig|2j and the self- 
trapping regime, i.e. z > Zc. Clearly the experimental 
setup allows the preparation of both dynamical regimes. 
After the initial preparation of the BEC in an asym- 
metric double-well potential, i.e. z y^ the dynamics is 
initiated by shifting the crossed dipole trap beam back 
to Ax = 0. The time scale of 7ms is non-adiabatic with 
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respect to the typical tunnelling time which is on the 
order of 50ms. In order to understand the experimental 
observations quantitatively this finite response time has 
to be taken into account (see 6 ). 



by analyzing the absorption image as shown in fig. [31 
From that measurement we deduce a lattice constant of 
di = 5.18(9)/im. 




Fig. 2 Preparation of the initial population imbalance. The 
filled red circles indicated the measured initial population im- 
balance corresponding to the two dynamical situations and 
the symmetric case. The condensate in the double-well po- 
tential is prepared by adiabatically ramping up the periodic 
potential. Since the symmetry of the effective double- well po- 
tential is given by the relative position Ax of the minimum of 
the crossed beam dipole trap and a maximum of the periodic 
potential the initial population can be adjusted by changing 
Ax. The graph reveals that the initial population imbalance 
can be created in the Josephson oscillation regime (indicated 
by the shading) and the self-trapping regime. The solid line 
is the theoretically predicted imbalance obtained by solving 
the non-polymomial Schrodinger equation numerically. 



The thorough calibration of all relevant parameters is 
the prerequisite for quantitative experiments. The har- 
monic confinement and the periodic potential are cal- 
ibrated independently. The harmonic trapping frequen- 
cies are deduced from standard frequency measurements 
of dipolar oscillations. They are excited by reducing the 
intensities of the laser beams within less than 500/is lead- 
ing to dipolar oscillations with small amplitudes (few 
jam) . These measurements lead to harmonic trapping fre- 
quencies of cOx = '^TT X 78(l)Hz, cOy ^ 27r X 90(l)Hz and 
cj^ ?^ 27r X 66(l)Hz. 

Since the tunnelling time is critically dependent on 
the thickness and height of the double- well barrier it is 
crucial to measure these parameters very accurately. As 
the width of the barrier is directly connected to the lat- 
tice spacing it can be measured in a straight forward 
way. A BEC of 2 x 10^ atoms is loaded into a deep opti- 
cal lattice {Vq ^ hx2 kHz) superimposed on a harmonic 
trap u;x,y,z = 27r x (5, 120, 140)Hz leading to a periodic 
array of pancake condensates (see fig. OJ. For these pa- 
rameters the influence of the shallow harmonic potential 
in x-direction on the spacing of the potential is negli- 
gible and thus the periodicity can be directly measured 
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Fig. 3 Calibration of the lattice spacing of the periodic po- 
tential, a) The absorption image of the atomic cloud reveals 
the periodic arrangement of pancake BECs. b) The den- 
sity profile in x-direction directly shows the lattice spacing 
di = 5.18(9)/im. From this measurement we can also deduce 
the optical resolution (sparrow-criterion) of our imaging sys- 
tem to be 3.2(2)/im. 



The potential height of the periodic potential is di- 
rectly connected to the barrier height. Due to the large 
well spacings the standard potential height calibration 
techniques do not work and thus we have developed a 
new method. The potential height is measured by ob- 
serving the relative motion of two BECs in the double- 
well potential. This motion can be excited by switching 
off the harmonic confinement in x-direction and ramping 
up the barrier height, i.e. the periodic potential height, 
by a factor of 5 within 2ms. This procedure leads to a 
non-adiabatic increase of the double- well potential spac- 
ing ddw = 4.2/im to the lattice constant of the peri- 
odic potential di = 5.18/im. The time scale of 2ms is 
chosen to reduce the collective excitations in the trans- 
verse direction but leads to very small oscillation ampli- 
tudes of approximately 400nm (see fig.EJ- However, this 
change of the center of mass separation can still be mea- 
sured with our optical imaging system which has an op- 
tical resolution of 3.2(2) /am as defined by the Sparrow- 
criterion. The experimental data are fitted with the re- 
sult obtained from the integration of the non-polynomial 
Schrodinger equation ^2) and a potential height of Vq = 
h X 412(20)Hz is deduced. This corresponds to a barrier 
height ofVb = hx 263 (20) Hz for the experimentally used 
harmonic confinement of cOx = ^tt x 78Hz. 



3 Josephson oscillations and self-trapping 

In contrast to all hitherto realized Josephson junctions 
in superconductors and superfluids, in our experiment 
the interaction between the tunnelling particles plays 
a crucial role. This nonlinear ity gives rise to new dy- 
namical regimes. Josephson oscillations, i.e. oscillation 
of population imbalance and relative phase of the two 
condensates, are predicted 13714,15 , if the initial pop- 
ulation imbalance of the two wells is below a critical 
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value. The dynamics changes drastically for initial pop- 
ulation differences above the threshold for macroscopic 
quantum self-trapping Zc n"6nT7il8 where large ampli- 
tude Josephson oscillations are inhibited and the phase 
difference increases in time. 

The experimental protocol for investigating the dy- 
namics is as follows: Rubidium atoms are precooled in a 
standard TOP trap, transferred into the crossed beam 
dipole trap and evaporatively cooled to degeneracy by 
lowering the light intensity. Finally the dipole laser beams 
and the optical lattice beams are ramped to the desired 
values. This sequence creates two weakly linked BECs 
inside an asymmetric double- well potential with well de- 
fined asymmetry and barrier height. The dynamics is 
initiated by shifting the crossed dipole beam realizing 
a symmetric double-well potential. After a given evolu- 
tion time the potential barrier is suddenly ramped up 
and the dipole trap beam is switched off. This results in 
dipole oscillations of the atomic cloud around two neigh- 
boring minima of the periodic potential as used for the 
calibration of the periodic potential height (see fig^. 
The atomic distribution is imaged at the time of max- 
imum separation using absorption imaging techniques. 
This protocol has been used for the first demonstra- 
tion of the transition between Josephson oscillations and 
macroscopic self-trapping in a single weak link 6 . In 
the following we shall show that the recent work by D. 
Ananikian and T. Bergeman 7 allows to understand the 
experimental observation quantitatively with a relatively 
simple two-mode model. 

The standard two-mode approximation assumes a 
weak link, i.e. the wave function overlap is negligible 
[TH| . For our experimental parameters this is not strictly 
true. S. Giovanazzi et al. ^Hj have already included the 
leading term of the correction to the simple two-mode 
model but recently D. Ananikian and T. Bergeman man- 
aged to include all terms and present their improved 
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Fig. 4 Calibration of the periodic potential height, a) The 
BEG is first loaded into the effective double-well potential 
(black solid line). Subsequently the harmonic confinement in 
x-direction is switched off and the periodic potential height is 
increased to Vb- This leads to dipole oscillations in the indi- 
vidual wells as indicated, b) The relative position of the two 
BEGs in the optical lattice as a function of time after ex- 
citation reveal the potential height. The experimental data 
(crosses) are compared to the numerical simulation (solid 
line), which has only the lattice height as a free parameter. 
We find a potential height of Vb = /i x 412(20)Hz. 



two- mode model in [7). They have derived differential 
equations for the basic two-mode parameters the pop- 
ulation difference z and the relative phase (\) between 
the two condensates. We will not elaborate on the the- 
oretical description any further and refer the reader to 
the reference 7 for details. Here, we report on the very 
good agreement of the prediction of this model with our 
experimental observation concerning the critical popu- 
lation imbalance distinguishing between the two differ- 
ent dynamical regimes. The simple constant tunnelling 
model predicts for our experimental parameters a critical 
population difference Zc = 0.23, but for this imbalance 
experimentally still Josephson oscillations are observed. 
From the experimental data we deduce z^^xp = 0.38(8) 
which is in very good agreement with the prediction of 
the improved two- mode model giving Zc = 0.35. This 
agreement is clearly revealed in fig|Sl where the phase 
plane portrait predicted by the time varying tunnelling 
model (solid lines) and the experimental data (circles) 
are shown. It is important to note that there are no free 
parameters in this graph, which is in contrast to the ear- 
lier reported phase-plane portrait 0, where the critical 
population difference was a free fit parameter. 
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Fig. 5 Gomparison of the experimentally obtained phase 
plane trajectories to the predictions of the extended two- 
mode model (solid lines). The Josephson oscillation regime 
(gray shaded region) is characterized by closed trajectories 
(filled circles). The separatrix, which is represented by the 
dashed line, constitutes the transition to the self-trapped 
regime (open circles). It is important to note that there are 
no free parameters. Thus the improved two-mode tunnelling 
model explains our observation quantitatively. 



Thus the observed dynamics can be understood in 
a simple two-mode model, but a time dependent tun- 
nelling rate has to be taken into account. This finding 
is an important prerequisite for further investigations of 
more complex systems which are build up by Josephson 
junctions. One interesting direction is the controlled re- 
alization of many coupled weak links in two dimensions 
where the topology could have a big influence on the 
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dynamics ^^. Also a very intriguing route is the inves- 
tigation of the influence of the residual thermal cloud 
on the coherent dynamics [21]. Here, a completely new 
regime can be reached with the BEC system since the 
thermal cloud and thus the dissipation can be controlled. 
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